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PLASTIC BUCKLING OF AXIALLY COMPRESSED 
ECCENTRICALLY STIFFENED CYLINDERS 

By Robert M. Jones 
Langley Research Center 

SUMMARY 

A solution for the plastic buckling of axially compressed eccentrically stiffened 
cylinders is derived by use of the J 2 deformation theory of plasticity (maximum energy 
of distortion theory) for a set of simply supported edge boundary conditions. Many effects 
which have not been studied previously are included, and a closed-form stability criterion 
is obtained. The present results are more realistic for practical stiffened cylinders than 
an existing solution. A numerical example is given to illustrate the effect of plastic action 
on the buckling load of a typical aerospace structure. 

INTRODUCTION 

The ring- and- stringer- stiffened cylinder is a commonly used aerospace structure 
for supporting axially compressive loads. Although considerable information is available, 
for determination of the elastic buckling load of such structures, information is scarce on 
buckling at stresses above the proportional limit of the shell material. Prior work on 
this problem is apparently confined to reference 1 wherein conventional orthotropic shell 
theory is employed to represent a stiffened cylinder. Two deficiencies occur in the 
approach of reference 1 ; the essentially uniaxial character of the stiffeners and the eccen- 
tricity of the stiffeners are not considered. 

In the present report, plastic buckling of axially compressed, simply supported, 
eccentrically stiffened cylinders such as shown in figure 1 is treated. The uniaxial char- 
acter of the stiffeners and the stiffener eccentricity (asymmetry about the cylinder middle 
surface) are explicitly accounted for by following the approach of references 2 and 3. 
Simple J 2 deformation theory of plasticity is utilized (also called maximum octahedral 
shear stress theory, maximum energy of distortion theory, etc.). Experimental results 
for buckling are in reasonably good agreement with Jg deformation theory. 

The present solution is an extension of Stowell's work (ref. 4) in which it is assumed 
that unloading does not occur during buckling. In contrast to reference 4, however, a 
variable Poisson's ratio is included in the present theory as was done by Bijlaard (ref. 5). 


Donnell-type stability differential equations are employed in the study, and it is assumed 
that buckling is a bifurcation from a membrane prebuckled shape. In addition, the cir- 
cumferential stiffeners are considered not loaded, that is, ring constraint is neglected, 
prior to buckling. An equation is developed which gives the buckling stress in terms of 
the geometry of the cylinder and the properties of the cylinder material. A numerical 
example is given to illustrate the use of the present solution and to show some typical 
results. 


SYMBOLS 


The units for the physical quantities used in this paper are given both in U.S. 
Customary Units and in the International System of Units (SI) (ref. 6). Factors relating 
these two systems are given in the appendix. 

a ring spacing (see fig. 1) 


A ij 

A r 


defined by equations (34) 


plasticity coefficients defined by equations (19) 


cross-sectional area of a ring 


As 

b 


cross-sectional area of a stringer 


stringer spacing (see fig. 1) 


B 


extensional stiffness, 


^sec ' 
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E sec 
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bending stiffness, 
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Young's modulus 
secant modulus 
tangent modulus 
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M x> M y } 

MxyjMy X 


| N x ,N y n 

N xy / 

N x 

R 


strain intensity (see eq. (2)) 

E 

shearing modulus, 

2(1 + Me) 

E 

secant shearing modulus, — sec 

2(1 + m) 

moment of inertia of stringer about middle surface of cylinder 

moment of inertia of ring about ring centroid 

moment of inertia of stringer about stringer centroid 

torsional constant of a ring 

torsional constant of a stringer 

material stress-strain curve parameter 

function defined by equation (14) 

cylinder length 

number of longitudinal buckle half-waves 
moments per unit length 

number of circumferential buckle waves 
in-plane forces per unit length 

applied axial force per unit length at buckling 
radius of cylinder middle surface 
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t 

u,v,w 

x,y,z 

z r 


cylinder thickness • 

displacements from a membrane prebuckled shape in longitudinal, cir- 
cumferential, and radial directions, respectively 

longitudinal, circumferential, and radial coordinates, respectively, on 
cylinder middle surface 

distance from ring centroid to cylinder middle surface (see fig. 1), 
positive when ring on outside 

distance from stringer centroid to cylinder middle surface (see fig. 1), 
positive when stringer on outside 

uniaxial strain 


e l> e 2» e 3 


variations in middle surface strain 


e x> € y>> / xy 


strains 


1 - 


E 


X = 


tan 
' sec 


4K 4 (1 - m 2 ) 




1 /l \ E 

Poisson's ratio, i - ju e ) — 


sec 


Me 


a 


a 


CT i 


V 

a oy 


Poisson's ratio in elastic state 

uniaxial stress 

axial buckling stress 

stress intensity defined by equation (1) 

stress at proportional limit 

0.2-percent offset yield stress 
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*■ 

^x^yj^xy stresses 

XijXojXo variations in middle surface curvature 

1 A O 

The prefix 6 denotes variation of the principal symbol during buckling. When 
the subscripts x or y follow a comma, they denote partial differentiation of the prin- 
cipal symbol with respect to x or y. Otherwise, the subscripts x and y denote 
the direction with which the principal symbol is associated. 

DERIVATION OF THEORY 


By use of simple J 2 deformation theory of plasticity, expressions are obtained 
for the variations of stresses during buckling in terms of the variations of strains and 
material properties during buckling. The reference surface is taken as the cylinder 
middle surface. The results are specialized for axially compressive loading. Subse- 
quently, the variations in stresses are integrated over the shell and stiffeners in order 
to obtain expressions for the variations in forces and moments during buckling. Finally, 
the variations in forces and moments are substituted in Donnell-type stability differential 
equations which are then solved to yield a closed-form stability criterion in terms of the 
geometric and material properties of the stiffened cylinder. 


Fundamental Relations in J 2 Deformation Theory 

The fundamental relations in the J 2 deformation theory of plasticity (see ref. 7) 
are the stress intensity = 3 t 0 /|/2, where t q is the octahedral shear stress): 


( 1 ) 


<*i = {°x 2 + CT y 2 “ CT x°y + 3rxy 2 ) 1//2 

the strain intensity e^ej = 3y 0 ^{2^2(1 + ju)J> where y Q is the octahedral shear strain j: 

1 /2 

ei = - - 2 j(l - M + M 2 )(ex 2 + e y 2 ) - (1 - 4/x + Ax* y + |U - # i ) 2 > / xy 2 j ( 2 ) 


and the stress-strain relations which are compatible with the stress and strain 
intensities: 
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° x = 7T3( 6x + “ ey ) 


E 


0y= fr^( ey + ,X£x )^ 


( 3 ) 


Esec 

xy = 2(T77o rxy 


In addition, certain material properties are defined in terms of the stress and strain 
intensities: first, the secant modulus E se c 


E 


7 i 


sec - g? 


(4) 


and second, the tangent modulus E tan 


®tan - 


dCTj 

dei 


(5) 


Poisson's ratio is a variable which is defined on page 387 of reference 7 as 

i E< 




'sec 


( 6 ) 


where /x e is the elastic value of /x. 


Variations of Strains and Stresses During Buckling 

During buckling, the strains vary from their prebuckling values. Let the variation 
be denoted by 6; then, 


6e x = e l " z Xx 
6e y = e 2 - zx 2 \ 

6 Xxy = 2e 3 " 2z *3 
J 


( 7 ) 


where €j, e 2 , and e 3 are the variations in middle surface strain and Xj> x 2 > and 
X 3 are the variations in middle surface curvature. 
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Corresponding variations in stresses can be obtained from equations (3); for 
example, 


6a x = + M5e y ) + ^e x + + e x E sec 6^— ^ 


M 


2 y 


Esec 5 -L (8) 


+ e v lL sec 


V 1 - 


Accordingly, the following variations are required: 

SEgec = - ^-fl - 


^i 


E 


sec/ 


(9) 


/ 1 \. 

2/x _ 

1 

n 

i 

i 1 - ^ 2 >1 


6/i 


(10) 



1 + M 2 

\1 - M 2 ) 

U - m 2 ) 2 _ 


6/1 


( 11 ) 


From equations (6) and (9), 


Hi -ft-feh 


( 12 ) 


so that all variations of the material parameters are expressed in terms of 6e^ which, 
in turn, can be obtained from the expression for (eq. (2)): 

5e\ = - 2(1 - /i + M 2 )( £ x6€x + £y6e y ) - (1 - 4 jli + /i 2 )(ex 5ey + 6y6€x) 

2K M ei(l - m 2)2L 

+ I<l - M) 2 rxy6yx^J (13) 
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where 


2 " Me E 


1 - u 2 E 


sec|i _ E tan 


+ 2(2 + fJ.)(J^Oy - 6(1 + 


Upon substitution of the stress-strain relations (eqs. (3)) and the variations of strains 
(eqs. (7)), equation (13) becomes 


J i rj— f (T2 - n)o x - (1 - 2/i)oJjej + [(2 - fi)cry - (1 - 2p)a^Je2 + 6(1 - M) T xy e 3 

2Kju° r i (l - M z ) \ 

- M) CT x - (1 - 2fi)a£j Xl + [(2 - M)a y - (1 - 2 m)ct x ]x 2 + 6(1 - MKxyXijjj ( 15 ) 


Finally, the relation for 6a x (eq. (8)) becomes 


60 v = 


— ^2 fl + /ie 2 " z (*l + ^2) 


E tan\ 1 




CT x + 9< ff y + 

1 - U* 


E sec 


x ^2 - p) CT x _ (1 “ 2fi)^£]ej + Q2 - M)cy - (1 - 2p.)a^je2 + 6(1 - p.) T xy e 3 

-{& - p)a x - (1 - 2/i)cry]x 1 + [(2 - p)cr y - (1 - 2 m)ctJx 2 + 6 (1 - ^TxyX^j (*6) 


The expression for 6a y is obtained by permutation of the x and y subscripts in 
equation (16) whereas the expression for Sixy is 
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6T xy - 


E 


sec 


2(1 - M 2 ) 


2(1 ■ M) ( £ S • **3) - + TT 


9 ■ Me E 


sec 


M E 


* |[(2 - n) cr x - (1 - 2pi)a^]e j + (J2 - n)(Jy - (1 - 2 p,)cr ^]€2 + 6(1 - p)t x y e 3 
- ^{[(2 ~ ^)°x ~ (1 “ 2fx)a^]xj + [(2 - - (1 - 2 /x)ct^]x 2 + 6(1 - M) T xyX3ji 


( 17 ) 


Equations (16) and (17) give the variations in stresses in terms of variations in cylinder 
middle surface strains and curvatures, membrane prebuckling stresses, and material 
properties in the prebuckled state. 

Hereafter, the equations are restricted to the case of axially compressive loading 
wherein Oj = a x because Oy = Txy = 0. The variations in stresses can now be written 
as 


*x = " Z *l) + ^ Al2 ( € 2 - 

6<Ty = ^^[ tAl2 ( ei ‘ ZX l) + A22 ( 6 2 ' Zx 2)]/ 

6T ^ = i(r^ A33 ( e3 ' ZX3 ) 

J 

where 

An = 1 - A(2 - fi) 2 
Al2 = 1 + - 2^)(2 - M) 

A22 = 1 - A(1 - 2fi) 2 
A33 = 2(1 - n) 


(18) 


(19) 
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in which 


1 - 


X = 


tan 


sec 


4K m (1 - m 2) 


( 20 ) 


For axial compression, reduces to 


k m = 


(i - |M<1 - 2rf 2 

1 | \ E sec/ 


4(1 - ju 2 ) 


(21) 


Variations of Forces During Buckling 

The variation of the forces during buckling is obtained by integration of the varia- 
tion of stresses over the appropriate area, for example, 


P */ 2 1 P 

6Ny = \ 60 v ^ + r \ dAs 

-t/2 b J 


-t/2 


( 22 ) 


which, upon substitution for 6 o x and integration, becomes 


6N x = rT^( Alle ! + ^ A 12e 2 ) + E tan - z sXl 


(23) 


The last term on the right-hand side of equation (23) is a result of the assumption that 
the stringers carry only axial stress (hence, the tangent modulus is appropriate). The 
expressions for 6Ny and 6Nxy are obtained in a similar manner; however, because 
the rings are assumed to carry no load (that is, ring constraint is neglected) prior to 
buckling, all ring terms have elastic moduli. In summary, 


-'V 


- E tan A s 


6N X = |BA n + E tan ~J e l + ^ BA 12 e 2 ' z s 

( EA r \ EA r 

6Ny = + BA 22 + —L U 2 - z r —L X2 




6N 


xy = (1 - m) B eg 


(24) 


where 


B = 


E 


sec 


1 - ju2 


(25) 
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• Variations of Moments During Buckling 

In a manner similar to that used to obtain the variation of the forces during 
buckling, the variation of the moments is obtained as 


SMv = - 


6M y = - 


6Mxy - 


i DAn + z^Etan + Etan -jHXi + M DA 12 X 2 " z s E tan e i 


( . — 2 EA r EI r \ _ EAj* 

mDA 12 X! + DA 22 + z r Z — + — )X 2 - z r — e 2 


(1 - m) D + Gsec T" 
o 


*3 


5Myx = — 


(l-M)D + 5^ 


*3 


} (26) 


where 


D = 


Esec t" 5 
12(1 - M 2 ) 


(27) 


Equations (24) and (26) for the variations in forces and moments during buckling reduce 
for the elastic eccentrically stiffened cylinder to equation (27) of reference 3. At this 
stage in the derivation, it is apparent that equations (11) and (12) of reference 1 apply 
only to orthotropic shells and, if applied to stiffened cylinders, imply that stiffeners have 
a two-dimensional stress state. The present theory, as does that of reference 3, treats 
the stiffeners as one -dimensional elements, a representation which is more realistic than 
an orthotropic shell model for practical stiffened cylinders. 


Stability Differential Equations 

The Donnell-type stability differential equations for axially compressed cylinders 

are 


-n 


*vr , __ n 

^x,x T wW xy,y u 


6N xy,x + 6Ny,y = 0 


■6M- 


x,xx + 6M xy,xy - 6%x,xy - 6 M y,yy + + N x w >xx 


8N V = 

+ «L + N V W YY = 0 




(28) 
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and the variations of middle surface strains and curvatures are 




e l = U ,X 

w 

Xj = w ,xx 

e 2 = v >y 

+ R 

*2 = w ,yy 

g — IjlZ. 

+ v 

X 3 = w ,xy^ 

€ 3 

2 


( 29 ) 


Upon substitution of the expressions for the variations in forces and moments during 
buckling from equations (24) and (26), respectively, and the variations of the middle sur- 
face strains and curvatures from equations (29), the stability differential equations 
(eqs. (28)) become 


BAn + Etan -j^j u ,xx + M&A.12 fv^yx + “ z sf E tan “j“^ w ,xxx + * ^ ^ ®( u ,yy + v ,xj) - 0 


^ 2 ^ B ( u ,y* + v > xx ) + ^ BA 12 u ,xy + (bA 2 2 + ^(v,yy + ^ - z r ^ w >y y 
^DAn + z s 2 Etan "j~^ + E tan r j~|j w ,xxxx + J^D^l " ^ + MA 12 ) + G sec 


= 0 


,xxyy 


+ IDA22 + 


Zj* EA r EI r 


_ EAr 


w 


a + ir/ w , yyyy ' z s E tan -£■ u ,xxx - z r -j^kyyy + “jpj 


/iBAjou x ( EA r 

+ ^^+(bA 2 2 + 




z r ~ a^” w ,yy + N x w ,xx = 0 


J 

(30) 


Stability Criterion 

It is desired to find the solution to the stability differential equations for the simply 
supported edge boundary conditions 

6N X = v = w = 6M X = 0 (31) 
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. The following buckling displacements satisfy the boundary conditions (eq. (31)): 


u = ucosH^cos^ 
l R 

v = 7sin“H S in 2Z 
l R 

w = wsinH^cos^ 
l R 


(32) 


(where u, v, and w are the amplitudes of the buckling displacements) and are sub- 
stituted in the stability differential equations (eqs. (30)). In order to obtain a nontrivial 
solution to the resulting equations, the determinant of the coefficients of u, v, and w 
must be zero, and the following stability criterion results: 



where 


•\ 




a 12 = 


nur n_ 
l R 


MBA12 + (^2^j ] 

f + ^ ta ^(x) 3 


mba 12 mjr 

a 13 = — 


a 22= M)B^r + (BA22 + 


a 23 = l BA 22 




EaA r EA r / n \ 3 

+ ~z—hr + z r — r— | =r I 


JR 


a IR 


a 33 = (^ DA 11 + z s 2 E t an + E tan + j * 1 * 1 - M + M a i 2 ) + G se c ^ ^ 

+ ( DA 22 + ^ 4 )( f ^ 4 ) 


(34) 
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The buckling stress is 


a = 


N x 



(35) 


Equation (33) is solved by trial and error for the plastic buckling load. Because of the 
numerous parameters in equation (33) and the need to investigate a large range of 
buckling modes (values of m and n), it is necessary from a practical standpoint to use 
a digital computer for numerical work. 

The solution represented by equation (33) reduces to: the solution of reference 8 
for unstiffened plastic cylinders, the classical Euler load for unstiffened elastic cylin- 
ders, and the solution of reference 2 for stiffened elastic cylinders. For a plastic sym- 
metrically stiffened plate column (z s = 0) with p = /i e = i, equation (33) reduces to 


N, 




^tan A os 


E sec ^ 


13 E tan | 
,4 4 E S ec/ 


whereas the results of reference 1 imply 


(36) 



Esec Iqs 
b 


3 Etanjos 

4 b 


+ 


Esec !\ 
9 \4 


+ 


3 Etan j 

4 E S ec/ 


(37) 


(See eqs. (53) and (54) of ref. 1.) Equation (37) contains a term associated with the stiff- 
eners which involves the secant modulus whereas intuition suggests that terms associated 
with stiffeners should involve only the tangent modulus in agreement with equation (36). 
Accordingly, equation (37) always gives higher plastic buckling stresses than does equa- 
tion (36). 


NUMERICAL EXAMPLE 

Because of the many geometric and material parameters in the theory, general 
results cannot be presented, even for a specific stress-strain relation. However, a 
numerical example is presented which is representative of current large -diameter 
booster interstage structures. 

The stress-strain curve for 7075-T6 aluminum alloy is shown in figure 2. This 
stress-strain curve was represented numerically by the Nadai stress-strain curve for a 
0.2-percent offset yield stress: 
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where for 7075-T6 aluminum alloy 

Upi = 55 ksi (0.38 GN/m 2 ) 

CT oy = 72 ksi (0.50 GN/m 2 ) 

E = 10.5 x 10 3 ksi (72 GN/m2) 

K = 2.2 

The cylinder has a radius of 135 in. (345 cm), a length of 135 in. (345 cm), and a thick- 
ness of 0.10 in. (0.25 cm). The stiffener dimensions are shown in figure 3; the stringer 
spacing is held constant at 2.00 in. (5.08 cm). 

For general instability (buckling in which both rings and stringers participate), 
buckling stresses are shown in figure 4 for a range of ring spacings. The increase in 
stress obtained by placing stringers externally as opposed to internally is seen to be con- 
siderably lessened for small ring spacings in this example. 

For panel instability (buckling in which only stringers participate, that is, buckling 
between rings), buckling stresses are shown in figure 5 for a range of ring spacings. 
Large effects of plastic action are seen in figure 5 for small ring spacings. However, 
only for large ring spacings are the panel instability stresses lower than the general 
instability stresses. The effects of eccentricity are seen to decrease with increasing 
plastic action. 


CONCLUDING REMARKS 

A solution for the plastic buckling of axially compressed eccentrically stiffened 
cylinders is derived by use of a deformation theory of plasticity for a set of simply sup- 
ported edge boundary conditions. Stiffeners are treated as one-dimensional elements, a 
representation which is more realistic than the orthotropic shell model employed in a 
previous plastic buckling analysis. A numerical example is given to illustrate the effect 
of plastic action on the buckling load of typical stiffened cylinders in aerospace 
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applications. For this example, the effects of eccentricity are seen to decrease with 
increasing plastic action. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., August 1, 1966, 

124-11-06-04-23. 


16 


APPENDIX 


CONVERSION OF U.S. CUSTOMARY UNITS TO SI UNITS 


The International System of Units (SI) was adopted by the Eleventh General 
Conference on Weights and Measures, Paris, October 1960, in Resolution No. 12 (ref. 6). 
Conversion factors for the units used herein are given in the following table: 


Physical 

U.S. Customary 

Conversion 

SI Unit 

quantity 

Unit 

factor 

(*) 

Length 

inch 

0.0254 

meters (m) 

Stress 

ksi = kipsf/in^ 

6.895 x 106 

newtons per square meter (N/m2) 


^Multiply value given in U.S. Customary Unit by conversion factor to obtain equiv- 
alent value in SI Unit. 

Prefixes to indicate multiple of units are as follows: 


Prefix 

Multiple 

centi (c) 
giga (G) 

10-2 

10 9 
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Figure 5.- Panel instability stresses as a function of a/R for axially compressed ring-and- stringer-stiffened cylinders. 
R = 135 inches (345 cm); t = 0.10 inch (0.25 cm); b = 2.00 inches (5.08 cm); 7075-T6 aluminum alloy. 
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